Influence of genetic variants and post transcriptional factors on imatinib transporters as determinants of the pharmacological response

Introduction – Although imatinib (IM) is a recognized gold standard in chronic myeloid leukemia (CML) therapy, resistance has emerged in a significant proportion of patients. Aim – The aim of this study was: (1) to investigate the role of genetic variants in genes encoding for IM transporters, as candidate of IM responsiveness and (2) to test the influence of miRNAs on IM response, focusing on efflux transporters. Methods – As a first step, a panel of polymorphisms (SNPs) was genotyped in a subgroup population of 189 patients enrolled in the Tyrosine Kinase Inhibitor Optimization and Selectivity (TOPS) trial. The association with cytogenetic response and molecular response (MR) was assessed for each SNP. As a second step, an in vitro IM-resistant model (K-562 CML cell line) was established. miRNAs profiles were analyzed using Taqman arrays and in silico search was performed for miRNAs deregulated after IM treatment. mRNA and protein expression were quantified using TaqMan realtime PCR and Western blotting, respectively. Results – (1) Among Caucasian patients, ABCB1 rs60023214 significantly correlated with complete MR (P = 0.005). Concerning SNPs combination in IM uptake transporters, the associations with treatment outcomes were statistically significant for both major and complete MR (P = 0.005 and P = 0.01, respectively). (2) ABCB1 protein was not expressed under any conditions of treatment, differently from ABCG2. Two deregulated miRNAs, namely miR-212 and miR-328, were identified to be inversely correlated with ABCG2 (r2= 0.57; p=0.03 and r2=0.47; p=0.06, respectively). Experiments of loss and gain of function confirmed the functional influence of these miRNAs on ABCG2. Conclusion – The multiple candidate gene approach identified single and combination of SNPs that can be proposed as predictor of IM response. The in vitro study suggested that IM resistance could be mediated by miRNA-dependent mechanism. Further studies are needed to validate these preliminary findings.

Pathogenic mechanisms in mitochondrial optic neuropathies

Leber’s hereditary optic neuropathy (LHON) and Autosomal Dominant Optic Atrophy (ADOA) are the two most common inherited optic neuropathies and both are the result of mitochondrial dysfunctions. Despite the primary mutations causing these disorders are different, being an mtDNA mutation in subunits of complex I in LHON and defects in the nuclear gene encoding the mitochondrial protein OPA1 in ADOA, both pathologies share some peculiar features, such a variable penetrance and tissue-specificity of the pathological processes. Probably, one of the most interesting and unclear aspect of LHON is the variable penetrance. This phenomenon is common in LHON families, most of them being homoplasmic mutant. Inter-family variability of penetrance may be caused by nuclear or mitochondrial ‘secondary’ genetic determinants or other predisposing triggering factors. We identified a compensatory mechanism in LHON patients, able to distinguish affected individuals from unaffected mutation carriers. In fact, carrier individuals resulted more efficient than affected subjects in increasing the mitochondrial biogenesis to compensate for the energetic defect. Thus, the activation of the mitochondrial biogenesis may be a crucial factor in modulating penetrance, determining the fate of subjects harbouring LHON mutations. Furthermore, mtDNA content can be used as a molecular biomarker which, for the first time, clearly differentiates LHON affected from LHON carrier individuals, providing a valid mechanism that may be exploited for development of therapeutic strategies. Although the mitochondrial biogenesis gained a relevant role in LHON pathogenesis, we failed to identify a genetic modifying factor for the variable penetrance in a set of candidate genes involved in the regulation of this process. A more systematic high-throughput approach will be necessary to select the genetic variants responsible for the different efficiency in activating mitochondrial biogenesis. A genetic modifying factor was instead identified in the MnSOD gene. The SNP Ala16Val in this gene seems to modulate LHON penetrance, since the Ala allele in this position significantly predisposes to be affected. Thus, we propose that high MnSOD activity in mitochondria of LHON subjects may produce an overload of H2O2 for the antioxidant machinery, leading to release from mitochondria of this radical and promoting a severe cell damage and death ADOA is due to mutation in the OPA1 gene in the large majority of cases. The causative nuclear defects in the remaining families with DOA have not been identified yet, but a small number of families have been mapped to other chromosomal loci (OPA3, OPA4, OPA5, OPA7, OPA8). Recently, a form of DOA and premature cataract (ADOAC) has been associated to pathogenic mutations of the OPA3 gene, encoding a mitochondrial protein. In the last year OPA3 has been investigated by two different groups, but a clear function for this protein and the pathogenic mechanism leading to ADOAC are still unclear. Our study on OPA3 provides new information about the pattern of expression of the two isoforms OPA3V1 and OPA3V2, and, moreover, suggests that OPA3 may have a different function in mitochondria from OPA1, the major site for ADOA mutations. In fact, based on our results, we propose that OPA3 is not involved in the mitochondrial fusion process, but, on the contrary, it may regulate mitochondrial fission. Furthermore, at difference from OPA1, we excluded a role for OPA3 in mtDNA maintenance and we failed to identify a direct interaction between OPA3 and OPA1. Considering the results from overexpression and silencing of OPA3, we can conclude that the overexpression has more drastic consequences on the cells than silencing, suggesting that OPA3 may cause optic atrophy via a gain-of-function mechanism. These data provide a new starting point for future investigations aimed at identifying the exact function of OPA3 and the pathogenic mechanism causing ADOAC.

Metodi computazionali per l'annotazione di genomi e proteomi

Il progresso tecnologico nel campo della biologia molecolare, pone la comunità scientifica di fronte all’esigenza di dare un’interpretazione all’enormità di sequenze biologiche che a mano a mano vanno a costituire le banche dati, siano esse proteine o acidi nucleici. In questo contesto la bioinformatica gioca un ruolo di primaria importanza. Un nuovo livello di possibilità conoscitive è stato introdotto con le tecnologie di Next Generation Sequencing (NGS), per mezzo delle quali è possibile ottenere interi genomi o trascrittomi in poco tempo e con bassi costi. Tra le applicazioni del NGS più rilevanti ci sono senza dubbio quelle oncologiche che prevedono la caratterizzazione genomica di tessuti tumorali e lo sviluppo di nuovi approcci diagnostici e terapeutici per il trattamento del cancro. Con l’analisi NGS è possibile individuare il set completo di variazioni che esistono nel genoma tumorale come varianti a singolo nucleotide, riarrangiamenti cromosomici, inserzioni e delezioni. Va però sottolineato che le variazioni trovate nei geni vanno in ultima battuta osservate dal punto di vista degli effetti a livello delle proteine in quanto esse sono le responsabili più dirette dei fenotipi alterati riscontrabili nella cellula tumorale. L’expertise bioinformatica va quindi collocata sia a livello dell’analisi del dato prodotto per mezzo di NGS ma anche nelle fasi successive ove è necessario effettuare l’annotazione dei geni contenuti nel genoma sequenziato e delle relative strutture proteiche che da esso sono espresse, o, come nel caso dello studio mutazionale, la valutazione dell’effetto della variazione genomica. È in questo contesto che si colloca il lavoro presentato: da un lato lo sviluppo di metodologie computazionali per l’annotazione di sequenze proteiche e dall’altro la messa a punto di una pipeline di analisi di dati prodotti con tecnologie NGS in applicazioni oncologiche avente come scopo finale quello della individuazione e caratterizzazione delle mutazioni genetiche tumorali a livello proteico.

Studio del possibile ruolo del complesso polycomb nel determinare la schisi del labbro e del palato

La labioschisi con o senza palatoschisi non-sindromica (NSCL/P) è tra le più frequenti alterazioni dello sviluppo embrionale, causata dall’interazione di fattori genetici e ambientali, moti dei quali ancora ignoti. L'obiettivo del mio progetto di Dottorato consiste nell’identificazione di fattori di rischio genetico in un processo a due stadi che prevede la selezione di geni candidati e la verifica del loro coinvolgimento nella determinazione della malformazione mediante studi di associazione. Ho analizzato alcuni polimorfismi a singolo nucleotide (SNPs) dei geni RFC1 e DHFR, appartenenti alla via metabolica dell’acido folico, evidenziando una debole associazione tra alcuni degli SNPs indagati e la NSCL/P nella popolazione italiana. Presso il laboratorio della Dott.ssa Mangold dell’Università di Bonn, ho valutato il ruolo di 15 diverse regioni cromosomiche nel determinare la suscettibilità alla malattia, evidenziando una significativa associazione per i marcatori localizzati in 8q24 e 1p22. Ho quindi rivolto la mia attenzione al ruolo del complesso Polycomb nell’insorgenza della schisi. Nell’uomo i due complessi Polycomb, PRC1 e PRC2, rimodellano la cromatina agendo da regolatori dei meccanismi trascrizionali alla base della differenziazione cellulare e dello sviluppo embrionale. Ho ipotizzato che mutazioni a carico di geni appartenenti a PRC2 possano essere considerati potenziali fattori di rischio genetico nel determinare la NSCL/P. Il razionale consiste nel fatto che JARID2, una proteina che interagisce con PRC2, è associata all’insorgenza della NSCL/P ed espressa a livello delle cellule epiteliali delle lamine palatine che si approssimano alla fusione. L’indagine condotta analizzando i geni di elementi o partner dei due complessi Polycomb, ha evidenziato un’associazione significativa con alcuni polimorfismi dei geni indagati, associazione ulteriormente confermata dall’analisi degli aplotipi. Le analisi condotte sui geni candidati mi hanno permesso di raccogliere dati interessanti sull’eziologia della malformazione. Studi indipendenti saranno necessari per poter validare l'associazione tra le varianti genetiche di questi geni candidati e la NSCL/P.

Studio sul ruolo dei fattori genetici coinvolti nella Valvola Aortica Bicuspide e/o nell'Aneurisma dell'Aorta Toracica

La Valvola Aortica Bicuspide (BAV) rappresenta la più comune anomalia cardiaca congenita, con un’incidenza dello 0,5%-2% nella popolazione generale. Si caratterizza per la presenza di due cuspidi valvolari anziché tre e comprende diverse forme. La BAV è frequentemente associata agli aneurismi dell’aorta toracica (TAA). La dilatazione dell’aorta espone al rischio di sviluppare le complicanze aortiche acute. Materiali e metodi Sono stati reclutati 20 probandi consecutivi sottoposti a chirurgia della valvola aortica e dell'aorta ascendente presso l'Unità di Chirurgia Cardiaca di Policlinico S.Orsola-Malpighi di TAA associata a BAV. Sono stati esclusi individui con una condizione sindromica predisponente l’aneurisma aortico. Ciascun familiare maggiorenne di primo grado è stato arruolato nello studio. L’analisi di mutazioni dell’intero gene ACTA2 è stata eseguita con la tecnica del “bidirectional direct sequencing”. Nelle forme familiari, l’intera porzione codificante del genoma è stata eseguita usando l’exome sequencing. Risultati Dopo il sequenziamento di tutti i 20 esoni e giunzioni di splicing di ACTA2 nei 20 probandi, non è stata individuata alcuna mutazione. Settantasette familiari di primo grado sono stati arruolati. Sono state identificate cinque forme familiari. In una famiglia è stata trovata una mutazione del gene MYH11 non ritenuta patogenetica. Conclusioni La mancanza di mutazioni, sia nelle forme sporadiche sia in quelle familiari, ci suggerisce che questo gene non è coinvolto nello sviluppo della BAV e TAA e, l’associazione che è stata riportata deve essere considerata occasionale. L’architettura genetica della BAV verosimilmente dovrebbe consistere in svariate differenti varianti genetiche che interagiscono in maniera additiva nel determinare un aumento del rischio.

Mechanism of resistance to tyrosine kinase inhibitors in philadelphia-positive acute lymphblastic leukaemia (all): from genetic alterations to impaired RNA editing

The Ph chromosome is the most frequent cytogenetic aberration associated with adult ALL and it represents the single most significant adverse prognostic marker. Despite imatinib has led to significant improvements in the treatment of patients with Ph+ ALL, in the majority of cases resistance developed quickly and disease progressed. Some mechanisms of resistance have been widely described but the full knowledge of contributing factors, driving both the disease and resistance, remains to be defined. The observation of rapid development of lymphoblastic leukemia in mice expressing altered Ikaros (Ik) isoforms represented the background of this study. Ikaros is a zinc finger transcription factor required for normal hemopoietic differentiation and proliferation, particularly in the lymphoid lineages. By means of alternative splicing, Ikaros encodes several proteins that differ in their abilities to bind to a consensus DNA-binding site. Shorter, DNA nonbinding isoforms exert a dominant negative effect, inhibiting the ability of longer heterodimer partners to bind DNA. The differential expression pattern of Ik isoforms in Ph+ ALL patients was analyzed in order to determine if molecular abnormalities involving the Ik gene could associate with resistance to imatinib and dasatinib. Bone marrow and peripheral blood samples from 46 adult patients (median age 55 yrs, 18-76) with Ph+ ALL at diagnosis and during treatment with imatinib (16 pts) or dasatinib (30 pts) were collected. We set up a fast, high-throughput method based on capillary electrophoresis technology to detect and quantify splice variants. 41% Ph+ ALL patients expressed high levels of the non DNA-binding dominant negative Ik6 isoform lacking critical N-terminal zinc-fingers which display abnormal subcellular compartmentalization pattern. Nuclear extracts from patients expressed Ik6 failed to bind DNA in mobility shift assay using a DNA probe containing an Ikaros-specific DNA binding sequence. In 59% Ph+ ALL patients there was the coexistence in the same PCR sample and at the same time of many splice variants corresponded to Ik1, Ik2, Ik4, Ik4A, Ik5A, Ik6, Ik6 and Ik8 isoforms. In these patients aberrant full-length Ikaros isoforms in Ph+ ALL characterized by a 60-bp insertion immediately downstream of exon 3 and a recurring 30-bp in-frame deletion at the end of exon 7 involving most frequently the Ik2, Ik4 isoforms were also identified. Both the insertion and deletion were due to the selection of alternative splice donor and acceptor sites. The molecular monitoring of minimal residual disease showed for the first time in vivo that the Ik6 expression strongly correlated with the BCR-ABL transcript levels suggesting that this alteration could depend on the Bcr-Abl activity. Patient-derived leukaemia cells expressed dominant-negative Ik6 at diagnosis and at the time of relapse, but never during remission. In order to mechanistically demonstrated whether in vitro the overexpression of Ik6 impairs the response to tyrosine kinase inhibitors (TKIs) and contributes to resistance, an imatinib-sensitive Ik6-negative Ph+ ALL cell line (SUP-B15) was transfected with the complete Ik6 DNA coding sequence. The expression of Ik6 strongly increased proliferation and inhibited apoptosis in TKI sensitive cells establishing a previously unknown link between specific molecular defects that involve the Ikaros gene and the resistance to TKIs in Ph+ ALL patients. Amplification and genomic sequence analysis of the exon splice junction regions showed the presence of 2 single nucleotide polymorphisms (SNPs): rs10251980 [A/G] in the exon2/3 splice junction and of rs10262731 [A/G] in the exon 7/8 splice junction in 50% and 36% of patients, respectively. A variant of the rs11329346 [-/C], in 16% of patients was also found. Other two different single nucleotide substitutions not recognized as SNP were observed. Some mutations were predicted by computational analyses (RESCUE approach) to alter cis-splicing elements. In conclusion, these findings demonstrated that the post-transcriptional regulation of alternative splicing of Ikaros gene is defective in the majority of Ph+ ALL patients treated with TKIs. The overexpression of Ik6 blocking B-cell differentiation could contribute to resistance opening a time frame, during which leukaemia cells acquire secondary transforming events that confer definitive resistance to imatinib and dasatinib.

TNFRSF13B Genetic variability an anthropological - evolutionary approach to Biomedical Research

In the recent years TNFRSF13B coding variants have been implicated by clinical genetics studies in Common Variable Immunodeficiency (CVID), the most common clinically relevant primary immunodeficiency in individuals of European ancestry, but their functional effects in relation to the development of the disease have not been entirely established. To examine the potential contribution of such variants to CVID, the more comprehensive perspective of an evolutionary approach was applied in this study, underling the belief that evolutionary genetics methods can play a role in dissecting the origin, causes and diffusion of human diseases, representing a powerful tool also in human health research. For this purpose, TNFRSF13B coding region was sequenced in 451 healthy individuals belonging to 26 worldwide populations, in addition to 96 control, 77 CVID and 38 Selective IgA Deficiency (IgAD) individuals from Italy, leading to the first achievement of a global picture of TNFRSF13B nucleotide diversity and haplotype structure and making suggestion of its evolutionary history possible. A slow rate of evolution, within our species and when compared to the chimpanzee, low levels of genetic diversity geographical structure and the absence of recent population specific selective pressures were observed for the examined genomic region, suggesting that geographical distribution of its variability is more plausibly related to its involvement also in innate immunity rather than in adaptive immunity only. This, together with the extremely subtle disease/healthy samples differences observed, suggests that CVID might be more likely related to still unknown environmental and genetic factors, rather than to the nature of TNFRSF13B variants only.

Analisi di “Copy Number Variants” ed identificazione di nuovi geni candidati per l’Autismo e Ritardo Mentale

I disturbi dello spettro autistico (DSA) ed il ritardo mentale (RM) sono caratterizzati da un’eziologia genetica complessa ed eterogenea. Grazie ai recenti sviluppi nella ricerca genomica, è stato possibile dimostrare il ruolo di numerose copy number variants (CNVs) nella patogenesi di questi disturbi, anche se nella maggior parte dei casi l’eziologia rimane ancora sconosciuta. Questo lavoro riguarda l’identificazione e la caratterizzazione dei CNVs in famiglie con DSA e RM. E’ stata studiata una microdelezione in 7q31 che coinvolge i geni IMMP2L e DOCK4, trasmessa dalla madre con dislessia a due figli con autismo ed una figlia con dislessia. Nella stessa famiglia segrega una seconda microdelezione in 2q14 che inattiva il gene CNTNAP5 ed è trasmessa dal padre (con tratti autistici) ai due figli con autismo. Abbiamo quindi ipotizzato che i geni DOCK4 e CNTNAP5 potessero essere implicati, rispettivamente, nella suscettibilità a dislessia e DSA. Lo screening di numerosi individui affetti ha supportato la nostra ipotesi, con l’identificazione di una nuova microdelezione di DOCK4 che segrega con la dislessia, e 3 nuove varianti missenso in CNTNAP5 in individui con autismo. Dall’analisi genomica comparativa su array (aCGH) di individui con RM, è stata identificata una delezione nella regione 7q31.32, che coinvolge il gene CADPS2, in due fratelli con RM e tratti autistici, probabilmente ereditata dalla madre. Lo screening di mutazione di questo gene in individui con autismo o RM, ha portato all’identificazione di 3 varianti non sinonime, assenti nei controlli, ed ereditate per via materna. Poiché CADPS2 risiede in una regione genomica che contiene loci soggetti ad imprinting, abbiamo ipotizzato che il gene CADPS2 possa essere anch’esso caratterizzato da imprinting, con espressione monoallelica materna. Lo studio di espressione di CADPS2 in cellule del sangue ha avvalorato questa ipotesi, implicando perciò CADPS2 come un nuovo gene di suscettibilità per il RM e DSA.

Investigating the role of Copy Number Variants in Specific Language Impairment and identification of new candidate genes

Specific language impairment (SLI) is a complex neurodevelopmental disorder defined as an unexpected failure to develop normal language abilities for no obvious reason. Copy number variants (CNVs) are an important source of variation in the susceptibility to neuropsychiatric disorders. Therefore, a CNV study within SLI families was performed to investigate the role of structural variants in SLI. Among the identified CNVs, we focused on CNVs on chromosome 15q11-q13, recurrently observed in neuropsychiatric conditions, and a homozygous exonic microdeletion in ZNF277. Since this microdeletion falls within the AUTS1 locus, a region linked to autism spectrum disorders (ASD), we investigated a potential role of ZNF277 in SLI and ASD. Frequency data and expression analysis of the ZNF277 microdeletion suggested that this variant may contribute to the risk of language impairments in a complex manner, that is independent of the autism risk previously described in this region. Moreover, we identified an affected individual with a dihydropyrimidine dehydrogenase (DPD) deficiency, caused by compound heterozygosity of two deleterious variants in the gene DPYD. Since DPYD represents a good candidate gene for both SLI and ASD, we investigated its involvement in the susceptibility to these two disorders, focusing on the splicing variant rs3918290, the most common mutation in the DPD deficiency. We observed a higher frequency of rs3918290 in SLI cases (1.2%), compared to controls (~0.6%), while no difference was observed in a large ASD cohort. DPYD mutation screening in 4 SLI and 7 ASD families carrying the splicing variant identified six known missense changes and a novel variant in the promoter region. These data suggest that the combined effect of the mutations identified in affected individuals may lead to an altered DPD activity and that rare variants in DPYD might contribute to a minority of cases, in conjunction with other genetic or non-genetic factors.

Analysis of Copy Number Variants identifies new candidate genes for Autism Spectrum Disorder and Intellectual Disability

Autism spectrum disorder (ASD) and Intellectual Disability (ID) are complex neuropsychiatric disorders characterized by extensive clinical and genetic heterogeneity and with overlapping risk factors. The aim of my project was to further investigate the role of Copy Numbers Variants (CNVs), identified through genome-wide studies performed by the Autism Geome Project (AGP) and the CHERISH consortium in large cohorts of ASD and ID cases, respectively. Specifically, I focused on four rare genic CNVs, selected on the basis of their impact on interesting ASD/ID candidate genes: a) a compound heterozygous deletion involving CTNNA3, predicted to cause the lack of functional protein; b) a 15q13.3 duplication containing CHRNA7; c) a 2q31.1 microdeletion encompassing KLHL23, SSB and METTL5; d) Lastly, I investigated the putative imprinting regulation of the CADPS2 gene, disrupted by a maternal deletion in two siblings with ASD and ID. This study provides further evidence for the role of CTNNA3, CHRNA7, KLHL23 and CADPS2 as ASD and/or ID susceptibility genes, and highlights that rare genetic variation contributes to disease risk in different ways: some rare mutations, such as those impacting CTNNA3, act in a recessive mode of inheritance, while other CNVs, such as those occurring in the 15q13.3 region, are implicated in multiple developmental and/or neurological disorders possibly interacting with other susceptibility variants elsewhere in the genome. On the other hand, the discovery of a tissue-specific monoallelic expression for the CADPS2 gene, implicates the involvement of epigenetic regulatory mechanisms as risk factors conferring susceptibility to ASD/ID.